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METHOD AND APPARATUS FOR ESTIMATING WATER BOTTOM 
REFLECTION TIME DIFFERENCES 



BACKGROUND OF THE INVENTION 
1. FIELD OF THE INVENTION 

This invention relates generally to marine seismic surveying, and, more particularly, 
to estimating water bottom time differences in marine seismic surveying. 



10 2. DESCRIPTION OF THE RELATED ART 

Seismic exploration is widely used to locate and/or survey subterranean geological 
formations for hydrocarbon deposits. Since many commercially valuable hydrocarbon 
deposits are located beneath bodies of water, various types of marine seismic surveys have 
been developed. In a typical marine seismic survey, such as the exemplary survey 

15 conceptually illustrated in Figure 1, an array 100 of marine seismic streamer cables 105 is 

towed behind a survey vessel 1 10 over a survey area 115. The seismic streamer cables 105 
may be several thousand meters long and contain a large number of sensors 125, such as 
hydrophones and associated electronic equipment, which are distributed along the length of 
the each seismic streamer cable 105. The survey vessel 110 also tows one or more seismic 

20 sources 120, such as airguns and the like. 

As the array 100 is towed over the survey area 115, acoustic signals, or "shots," 
produced by the seismic sources 120 are directed down through the water into the earth 
beneath (not shown), where they are reflected from the various subterranean geological 
25 formations. The reflected signals are received by the sensors 125 in the seismic streamer 

cables 105, digitized and then transmitted to the survey vessel 110. The digitized signals are 
referred to as "traces" and are recorded and at least partially processed at the survey vessel 
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110. The ultimate aim of this process is to build up a representation of the subterranean 
geological formations beneath the array 100. Analysis of the representation may indicate 
probable locations of hydrocarbon deposits in the subterranean geological formations. 

Since the area of the array 100 is typically much smaller than the survey area 115, a 
representation of the earth strata in the survey area 1 1 5 may be formed by combining data 
collected along a plurality of sail lines 130(l-n). For example, a single survey vessel 110 
may tow a single array 100 along each of the sail lines 130(l-n). Alternatively, a plurality of 
survey vessels 110 may tow a plurality of arrays 100 along a corresponding plurality of the 
sail lines 130(l-n). However, variations in the water conditions, e.g. water temperature, 
salinity, and the like, between the plurality of sail lines 130(l-n) may cause variations in the 
velocity of sound in water among the sail lines 130(l-n). For example, the variations in 
seismic travel time can be on the order of 10 or 20 milliseconds for traces having a small 
distance between the source and detector for surveys carried out in deeper waters (greater 
than 200 m). The variations in the seismic wave travel times may shift the temporal position 
of the various events recorded in the seismic data, including, but not limited to, reflections 
and refractions of the seismic waves from the subterranean geological formations beneath the 
array 100. Consequently, the variations in the travel times may make it difficult to analyze 
the combined seismic data set and may reduce the accuracy of the survey. 

Moreover, the data for the sail lines 130(l-n) may be collected at different times. For 
one example, a single pass along one of the sail lines 130(l-n) may take several hours to 
complete so, if a single survey vessel 1 10 is used, data for the first sail line 130(1) will be 
recorded at an earlier time than data for the last sail line 130(n). For another example, 
inclement weather and/or high seas may force a survey to be suspended before resuming 
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hours or days later. For yet another example, historical data from previous surveys performed 
months or years earlier may be combined with new data to extend the survey or to fill in 
deficiencies in coverage that may be introduced by currents, obstacles such as platforms, and 
the like. And for yet another example, data from repeat surveys may be used to analyze and 
5 monitor changes in productive oil and/or gas reservoirs. 

Combining data from different times, and especially from different surveys, may 
exacerbate the aforementioned difficulties associated with variations in the velocity of sound 
in the water layer. For example, seasonal variations of the water temperature, salinity, and the 
10 like, may cause pronounced variations in the velocity of sound in water. For another example, 

shifts in water currents may cause unpredictable variations in the velocity of sound in water, 
particularly for surveys carried out near the edge of strong water currents. 

The seismic data is commonly corrected for the variations in the velocity of sound in 
15 water by computing one or more so-called "delta t" (At) values, which are typically defined 

as a difference between an expected travel time, usually based on an assumed ideal water 
velocity, and a measured travel time for one or more seismic signals. For example, the 
assumed ideal water velocity may be a constant velocity or one with very smooth spatial 
changes in velocity. 

20 

In one conventional method of determining the delta-t values (described in Wombell, 
R., 1996, "Water velocity variations in 3-D seismic processing," 66th Ann. Internat. Mtg: 
Society of Exploration Geophysicists, Expanded Abstracts, 1666-1669.), normal move-out 
stacking velocities and zero-offset water bottom times are computed along adjacent sail lines. 
25 The velocities are then converted to zero offset travel time differences using the formula: 
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At = T W (AV W /V W ), where At is the difference in two-way travel time at zero offset due to the 
change in water velocity, T w is the zero offset water bottom time, V w is the reference water 
velocity chosen by the practitioner, and AV W is the difference between V w and the computed 
stacking velocity. The delta-t values are then applied to normal move-out corrected seismic 
data. One problem with this method is that the velocity analysis must be extremely accurate. 

i 

Another issue is the effect of water bottom structure on the velocity analysis. If the dip of the 
water bottom (the angle the water bottom makes with a horizontal plane) changes between or 
along sail lines, the calculated velocities are strongly affected and may reduce the accuracy of 
the At calculation. 

Another method of determining the delta-t values that may be used to form combined 
data sets is described in Fried, J., and MacKay, S., 2001, "Dynamic Corrections for Water 
Velocity Variations: a Nova Scotia case history," Canadian Society of Exploration 
Geophysicists, October 2001 technical luncheon. In this method, normal move-out 
corrections are applied to pick times with a single velocity. The corrected pick times are then 
averaged for each combination of sail line, cross line, and common midpoint. The 
differences between sail line-cross line groups having overlapping midpoints are then 
evaluated and reduced using an iterative method. The method produces delta-t values that are 
used to apply a dynamic correction to the seismic data. One difficulty with this method is that 
the averaged pick times are affected by the difference between the actual normal move-out 
and the approximate normal move-out applied to correct the pick times. The effect of these 
differences propagates into the delta-t values. Reducing the offset range of pick times used in 
the average reduces the differences between the actual normal move-out and the approximate 
normal move-out applied to correct the pick times. However, reducing the offset range may 
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also reduce the amount of overlapping data upon which the method depends. Also, since the 
move-out is affected by the dip of the water bottom, changes in dip between sail lines can 
also affect the delta-t values. Furthermore, the iteration procedures used in this method are 
difficult to apply in practice. 

5 

SUMMARY OF THE INVENTION 

In one aspect of the instant invention, a method is provided for estimating water 
bottom reflection time differences. The method includes selecting a plurality of trace groups, 
each trace having an offset from a midpoint and a water bottom reflection time. The method 

10 also includes determining a rate of change of the water bottom reflection time of each trace 

and generating a model water velocity and a model water bottom reflection time for each 
trace midpoint using a preselected function of the water bottom reflection times and the rates 
of change. Furthermore, the method includes generating a plurality of water bottom reflection 
time corrections for the traces in each trace group using a pre-selected function of the water 

1 5 bottom reflection times, the model water bottom reflection times, the model water velocities, 

and the trace offsets. In other aspects of the present invention, an article comprising one or 
more machine-readable storage media may contain instructions that when executed enable a 
computer to perform the method or, alternatively, the article may contain data and/or data 
structures formed by the method. 

20 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention may be understood by reference to the following description taken in 
conjunction with the accompanying drawings, in which like reference numerals identify like 
elements, and in which: 

25 
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Figure 1 conceptually illustrates an exemplary prior art marine seismic survey; 

Figure 2 conceptually illustrates two sail lines of a marine seismic survey area; 

Figures 3A and 3B conceptually illustrate a survey vessel at two locations, as well as 
acoustic signals generated by the survey vessel, and reflected signals received by the survey 
vessel; 

Figure 4 conceptually illustrates a plurality of midpoint cells in the marine seismic 
survey area shown in Figure 2; 

Figure 5 conceptually illustrates a plurality of signal paths from a surface through a 
midpoint on a sea floor and back to the surface; 

Figure 6 conceptually illustrates a method for analyzing traces collected along a 
plurality of sail lines, such as the sail lines shown in Figure 2; and 

Figures 7 A and 7B conceptually illustrate aspects of a computing apparatus that may 
be used to implement the method shown in Figure 6. 

While the invention is susceptible to various modifications and alternative forms, 
specific embodiments thereof have been shown by way of example in the drawings and are 
herein described in detail. It should be understood, however, that the description herein of 
specific embodiments is not intended to limit the invention to the particular forms disclosed, 
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but on the contrary, the intention is to cover all modifications, equivalents, and alternatives 
falling within the spirit and scope of the invention as defined by the appended claims. 

DETAILED DESCRIPTION OF SPECIFIC EMBODIMENTS 

5 Illustrative embodiments of the invention are described below. In the interest of 

clarity, not all features of an actual implementation are described in this specification. It will 
of course be appreciated that in the development of any such actual embodiment, numerous 
implementation-specific decisions must be made to achieve the developers' specific goals, 
such as compliance with system-related and business-related constraints, which will vary 
10 from one implementation to another. Moreover, it will be appreciated that such a 

development effort might be complex and time-consuming, but would nevertheless be a 
routine undertaking for those of ordinary skill in the art having the benefit of this disclosure. 

Figure 2 conceptually illustrates a marine seismic survey area 200. To survey the 
15 marine seismic survey area 200, one or more survey vessels 210(1-2) tow one or more 

seismic arrays 215(1-2) over the marine seismic survey area 200. It will also be appreciated 
that, while the survey vessels 210(1-2) typically operate on the surface of the sea, the marine 
seismic survey area 200 refers to a portion of the sea bed. Furthermore, the present invention 
is not limited to undersea exploration, and may also be applied to surveys undertaken in 
20 freshwater, brackish water, and the like. 

In the illustrated embodiment, the seismic arrays 215(1-2) include a plurality of 
seismic sources 120 and seismic sensors 125, such as hydrophones, geophones, and the like, 
which may be coupled to the survey vessel 210(1-2) by cables 105. The plurality of seismic 
25 sources 120 and seismic sensors 125 may communicate with a signal processing unit 230 
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deployed on the survey vessels 210(1-2). In one embodiment, the plurality of seismic sources 
120 and seismic sensors 125 may communicate with the signal processing unit 230 via the 
cables 105. For example, the cables 105 may include wires, fiber-optic cables, and the like 
that may allow the seismic sources 120 and the seismic sensors 125 to signal processing unit 
230 to exchange signals. 

However, the seismic sensors 125 and seismic sources 120 may not always be 
t deployed as a part of the seismic arrays 215(1-2). In some alternative embodiments, the 
seismic sensors 125 can be deployed on the ocean bottom instead of being towed behind the 
survey vessels 210(1-2). For example, the seismic sensors 125 may be deployed on one or 
more ocean-bottom cables ("OBC"). The OBCs may be deployed on the seafloor so that the 
seismic sensors 125 may record and relay data to the signal processing units 230 on the 
seismic survey vessels 210(1-2). Furthermore, in one alternative embodiment, the seismic 
sources 120 may be deployed on buoys (not shown). In another alternative embodiment, the 
seismic sources 120 may be towed by a second vessel (not shown). 

The survey vessels 210(1-2) tow the seismic arrays along a plurality of sail lines, such 
as the two sail lines 220(1-2). Although only two sail lines 220(1-2) are shown in Figure 2 A, 
persons of ordinary skill in the art will appreciate that surveying the marine seismic survey 
area 200 typically requires more than two sail lines 220(1-2). For example, a survey covering 
an area of 40 x 70 miles requires about 160 sail lines 220(1-2), with each sail line 220(1-2) 
capturing about 1300 feet of subsurface coverage perpendicular to the direction of boat travel. 
Furthermore, although the two sail lines 220(1-2) shown in Figure 2 A are depicted as being 
approximately adjacent, those of ordinary skill in the art will appreciate that this is not 
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necessary for the practice of the present invention. In one embodiment, the two sail lines 
220(1-2) may not be adjacent or, alternatively, the two sail lines 220(1-2) may overlap. 

Figures 3 A and 3B conceptually illustrate a side view of the survey vessel 210(1) and 
a portion of the seismic array 215(1) at two different locations. In operation, the seismic 
source 125 shown in Figure 3 A provides an acoustic signal 300(1) that propagates to a 
seismic sensor 310(1) through a reflection point 320 located on a sea floor 325 between the 
seismic source 125 and the seismic sensor 310(1). Similarly, the seismic source 125 provides 
an acoustic signal 300(2) that propagates to a seismic sensor 310(2) through a reflection point 
330, as shown in Figure 3B. In the illustrated embodiment, the sea floor 325 is flat and so the 
reflection points 320, 330 are half-way between the seismic source 125 and the seismic 
sensors 310(1-2). However, those of ordinary skill in the art will appreciate that the reflection 
points 320, 330 may not necessarily be located half-way between the seismic source 125 and 
the seismic sensors 310(1-2). For example, a dipping sea floor 325 may change the location 
of the reflection points 320, 330. 

In one embodiment, signals are generated by the seismic sensors 310(1-2) in response 
to receiving the reflected and/or refracted acoustic signals 300(1-2) and then the generated 
signals are transmitted to the signal processing unit 230 on the survey vessel 210(1-2). In the 
illustrated embodiment, the signal processing unit 230 uses the transmitted signals to form 
one or more traces representative of the reflected and/or refracted acoustic signals 300(1-2), 
in a manner well known in the art. The signals and/or the traces may be stored on any storage 
medium, including, but not limited to, recording tape, magnetic disks, compact disks, and 
DVDs. Some embodiments may also, in addition to or in lieu of storing the signals and/or the 
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traces, transmit them to an on-shore facility. This may be done, for example, over a satellite 
link. 

The traces may be grouped according to the location of the reflection points 320, 330. 
5 For example, in the embodiment illustrated in Figure 4, the marine seismic survey area 200 is 

divided into a plurality of midpoint cells 401. The signals provided by the seismic array 
215(1) (not shown in Figure 4) during a single pass over the marine seismic survey area 200 
have reflection points 320 that are distributed in a band 405(1) of midpoint cells 401. 
Similarly, the signals provided by the seismic array 215(2) (not shown in Figure 4) during a 
10 single pass over the marine seismic survey area 200 may have reflection points 330 that are 

distributed in a band 405(2) of midpoint cells 401. Although the bands 415 (1-2) shown in 
Figure 4 are depicted as adjacent, this is not necessary for the practice of the present 
invention, as discussed above. 

15 Traces having a common midpoint cell 401 may be combined into a single data set, a 

process known in the art as bringing the traces to a common midpoint, forming a common 
midpoint gather, and the like. Furthermore, the midpoint cells 401 may be combined into 
cross-line groups 415(1-2). In one embodiment, traces corresponding to the midpoint cells 
401 in the cross-line group 415(1) may then be combined to form a first sail-line/cross-line 

20 gather. Similarly, traces corresponding to the midpoint cells 401 in the cross-line group 

415(2) may then be combined to form a second sail-line/cross-line gather. In one 
embodiment, traces associated with different sail lines 220(1-2) may also be combined. For 
example, the traces corresponding to the midpoint cells 401 in the cross-line groups 415(1-2) 
may all be combined to form a cross-line gather. 

25 
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The gathered traces may then be stacked to form a representation of various structures 
(not shown) within and/or above the seabed in the subsurface marine seismic survey area 
200, in a manner well known to those of ordinary skill in the art. However, as discussed 
above, variations in the water conditions, e.g. water temperature, salinity, and the like, may 
5 cause variations in the velocity of sound in water among the traces in the combined data sets. 

In accordance with common usage in the art, the velocity of sound in water will hereinafter 
be referred to as the "water speed" or the "water velocity." Variations in the water velocity 
among the traces in a combined data set, or gather, may make it difficult to analyze the 
combined data sets and may reduce the accuracy of the survey. 

10 

Figure 5 conceptually illustrates a plurality of signal paths 500(1-3) from a surface 
505 through a common midpoint cell 510 at a sea floor 520 and back to the surface 505. Each 
signal path 500(1-3) has a corresponding offset 530(1-3), which represents a horizontal 
separation of the seismic source and seismic sensor. Thus, for a water layer 540 having a 
1 5 depth Z w and a velocity of sound in water, or a water velocity, of V w , an acoustic signal that 

propagates along the signal paths 500(1-3) has a corresponding travel time Tj. 3 given by the 

formula T x _ 3 = (T 0 2 +X X _ 3 IVj) X12 , where Xjs is the length of the corresponding offset 
530(1-3) and T 0 = 2Z W /V W is the vertical two-way travel time, i.e. the travel time of an 
acoustic signal propagating along the line 550. The vertical two-way travel time is also 
20 referred to hereinafter as the zero-offset travel time. 

When combining traces formed from signals that propagate along the signal paths 
500(1-3), it is conventional to apply a so-called normal move-out (NMO) correction to the 
traces. The NMO correction includes transforming a time coordinate of the traces using the 
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equation r 0 I _ 3 = (T x _ 2 - X x _ 2 /V w 2 ) ]/2 . If the water velocity V w is the same for all the traces, 

then the NMOcorrected travel times To, 1-3 are all equal to a zero-offset travel time To. 
Grouping and/or combining the NMOcorrected traces often improves the representation of 
the marine seismic survey area 200 by, e.g., increasing the signal-to-noise ratio of the data. 
5 However, as discussed above, variations in the water conditions, e.g. water temperature, 
salinity, and the like, between the plurality of sail lines 200(1-2) may cause sail-line-to-sail- 
line variations in the water velocity, such that the NMOcorrected travel times Toj.3 are 
different for different traces. 

10 Figure 6 conceptually illustrates a method for analyzing traces collected along a 

plurality of sail lines (e.g., the sail lines 220(1-2) shown in Figure 2) in a marine seismic 
survey area 200. Groups of traces are selected (at 600). In one embodiment, selecting (at 600) 
the trace groups includes gathering traces having a common midpoint cell 410 in a common 
cross-line group 415(1-2), as discussed above. However, in alternative embodiment, any 

15 desirable grouping of the traces may be used. For example, traces having, or expected to 

have, a common water velocity may be grouped together. 

In the illustrated embodiment, the rate of change of the water bottom reflection time 
for each trace is determined (at 610). For example, the water bottom reflection time, T m may 
20 be determined by NMOcorrecting pick times, T, associated with each trace, according to the 

formula: 

(l) T n ={r 2 -x 2 s 2 )'\ 

where X is the offset of the trace and S is the inverse of the water velocity. The inverse of the 
water velocity is commonly referred to as the "slowness." The rate of change of the water 
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bottom reflection time is then determined (at 610) by taking the derivative of equation (1) 
with respect to the slowness. Thus, in the illustrated embodiment, the rate of change of the 
water bottom reflection time with respect to the slowness, P, is given by the formula: 

(2) j3 = ST„/SS = -X 2 S/T n 
However, it will be appreciated that the rate of change of the water bottom reflection time, J3, 
is given in equation (2) is merely intended to be illustrative and not to limit the present 
invention. In alternative embodiments, other expressions for the rate of change of the water 
bottom reflection time, including derivatives with respect to the water velocity and the like, 
may be used. 

Model water velocities and model water bottom reflection times are then generated (at 
620). In one embodiment, the rate of change of the water bottom reflection time may be used 
to define a function of the NMO-corrected travel times, T n (i,j), and a plurality of model 
travel times, T 0 @), associated with the midpoint cells, e.g., the midpoint cells 410 shown in 
Figure 4, For example, the function may be a square difference between the NMO-corrected 
travel times, T„(ij) and the plurality of model travel times, T 0 (j)\ 

(3) e 2 =(T 0 (j)+/3AS(j)-T„(i,j)Y 

In equation (3), the index "i" refers to the trace, the index "j" refers to the midpoint cell, and 
the rate of change of the water bottom reflection time with respect to the slowness, p, is given 
by equation (3). The model water velocity is included in equation (3) via the model 
slowness, AS(j) for the midpoint cell indicated by the index "j". 

In some alternative embodiments, a spatially smooth model water velocity is 
computed. For example, the model slowness, AS(j), may be an average taken over adjacent 
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midpoint cells. In one embodiment, the average may be a 3-cell average in which AS(j) is 
replaced with AS(jj), where the index jj indicates the group of three cells, -j-L, j, 
Alternatively, a weighted average based on the number of traces falling in cells j-l, j, j+1, 
could be taken. In yet another alternative embodiment, AS(j) may be averaged and/or 
smoothed after performing the minimization. 

The model water velocities and model water bottom reflection times may be 
generated (at 620) by minimizing the squared difference given by equation (3) using, e.g. a 
least-squares solution such as a Gauss-Seidel method or a conjugate gradient method. 
Although not necessary for the practice of the present invention, in various alternative 
embodiments, the above process may be iterated to generate iteratively updated model water 
velocities and iteratively updated model water bottom reflection times. 

The model water velocities extracted during the above-described velocity analysis 
may be affected by dip of the water bottom and the orientation of the line connecting the 
source and detector (or source-detector azimuth), as will be appreciated by those of ordinary 
skill in the art. Thus, in alternative embodiments, the pick times may be normalized using the 
azimuth and an approximate dip field measured from common midpoint stacks of the water 
bottom. However, normalizing the pick times is not necessary for the practice of the present 
invention and may be omitted. Moreover, those of ordinary skill in the art will appreciate 
that, in alternative embodiments, other selected portions of the seismic data and/or alternative 
techniques may also be used to generate (at 620) the model water velocity and/or the model 
water bottom reflection times. The alternative techniques include, but are not limited to, other 
averaging techniques, estimation procedures, direct measurements, and the like. These 
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techniques may be applied to the pick times and/or the other selected portions of the seismic 
data. 

One or more reflection times, such as the zero-offset reflection times, To(j) ? and one 
or more water bottom reflection time corrections, such as the zero-offset water bottom 
reflection time corrections, Ato(k), are then generated (at 630). In one embodiment, 
generating (at 630) the zero-offset reflection times and zero-offset water bottom reflection 
time corrections includes determining NMO-corrected travel times T n (iJ,k) by NMO- 
correcting the pick times T(i,j,k) using the model water velocities, represented here by the 
slowness S(j), determined in the aforementioned velocity analysis and the equation: 

(4) s T H (iJ 9 k) = (T(i,J,k) 2 -X(iJ,k) 2 S(j) 2 ) m 

The index "k" in equation (4) refers to the selected group. In one embodiment, the NMO- 
corrected travel times T n (iJ,k) are modeled as a linear combination of zero-offset travel times 
T 0 (j) and a plurality of zero-offset time corrections Ato(k) 9 as in the equation: 

(5) T n (/, j, k) = T 0 U) + a(i, h *) A/ 0 (*) 

where a(ij,k) is a function that relates zero-offset time corrections to non-zero offset time 
corrections. In one embodiment, a(i,j,k) 9 is given approximately by the equation: 

(6) a(h j 9 k) » 1 + (X(i 9 j, k)s(j)/ T„ (/, k)) 2 . 

However, it will be appreciated by those of ordinary skill in the art that the functional 
definition of a(i j,k) is a matter of design choice and not material to the present invention. 

Equation (5) may then be solved for the updated zero-offset travel times To(j) and a 
plurality of zero-offset time corrections Ato(kk) using a variety of techniques well known to 
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those of ordinary skill in the art. For example, equation (5) the solved by minimizing e 1 in 
the equation: 

(7) ^^(^(^y^J-roOO+aO^/.^AroWy- 

The least-squares error function defined in equation (7) may be solved for the zero- 
offset travel time corrections Ato(k) that minimize the error using a Gauss-Seidel method. 
However, those of ordinary skill in the art will appreciate that equations (7) may be solved by 
a variety of techniques. For example, in one alternative embodiment, equation (7) may be 
solved using a conjugate gradient method. In one other alternative embodiment, additional 
constraints may be added to equation (7) to facilitate a solution. In yet another alternative 
embodiment, iterations of the minimization of equation (7) may be performed using updated 
values of T n and a computed from the T n and a sum of the At 0 (k) of the previous "m" 
iterations: 

(8) T n (i, j, k 9 m + l) = T n (/, j, k,0) + X At 0 (*, m) 
until a selected convergence criterion is reached. For example, the minimizations may be 
iterated until the changes in the water bottom correction times At 0 (k,m) change in absolute 
value by less than 1 percent of the absolute value of the values from the first iteration, 
At 0 (k,0). 

The water bottom times and water bottom reflection time corrections (To(j) and 
At 0 (k), respectively) are then used to update (at 640) the water bottom reflection times. For 
example, in one embodiment, the water bottom reflection times may be updated (at 640) 
using a dynamic correction procedure, known to those of ordinary skill in the art. However, 
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in alternative embodiments, any desirable procedure for updating (at 640) the water bottom 
reflection times may be used. 

Figures 7A and 7B show a computing apparatus 700 that may be used to perform the 
aforementioned operations. The computing apparatus 700 includes a processor 705 
communicating with some storage 710 over a bus system 715. The storage 710 may include 
a hard disk and/or random access memory ("RAM") and/or removable storage such as a 
floppy magnetic disk 717 and an optical disk 720. The storage 710 is encoded with a data 
structure 725 storing the signals collected as discussed above, an operating system 730, user 
interface software 735, and an application 765. The user interface software 735, in 
conjunction with a display 740, implements a user interface 745. The user interface 745 may 
include peripheral I/O devices such as a key pad or keyboard 750, a mouse 755, or a joystick 
760. The processor 705 runs under the control of the operating system 730, which may be 
practically any operating system known to the art. The application 765 is invoked by the 
operating system 730 upon power up, reset, or both, depending on the implementation of the 
operating system 730. 

As discussed above, data collected during the marine seismic survey may be 
communicated to the computing apparatus 700 via any storage medium, including, but not 
limited to, recording tape, magnetic disks, compact disks, and DVDs. The data collected 
during the marine seismic survey may also be communicated directly to the computing 
apparatus 700 by, e.g., a satellite link 770, and stored in the storage 710. Some portions of 
the detailed descriptions herein are consequently presented in terms of a software 
implemented process involving symbolic representations of operations on data bits within a 
memory in a computing system or a computing device. These descriptions and 
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representations are the means used by those in the art to most effectively convey the 
substance of their work to others skilled in the art. The process and operation require 
physical manipulations of physical quantities. Usually, though not necessarily, these 
quantities take the form of electrical, magnetic, or optical signals capable of being stored, 
transferred, combined, compared, and otherwise manipulated. It has proven convenient at 
times, principally for reasons of common usage, to refer to these signals as bits, values, 
elements, symbols, characters, terms, numbers, or the like. 

It should be borne in mind, however, that all of these and similar terms are to be 
associated with the appropriate physical quantities and are merely convenient labels applied 
to these quantifies. Unless specifically stated or otherwise as may be apparent, throughout 
the present disclosure, these descriptions refer to the action and processes of an electronic 
device, that manipulates and transforms data represented as physical (electronic, magnetic, or 
optical) quantities within some electronic device's storage into other data similarly 
represented as physical quantities within the storage, or in transmission or display devices. 
Exemplary of the terms denoting such a description are, without limitation, the terms 
"processing," "computing," "calculating," "determining," "displaying," and the like. 

Note also that the software implemented aspects of the invention are typically 
encoded on some form of program storage medium or implemented over some type of 
transmission medium. The program storage medium may be magnetic (e.g., a floppy disk or 
a hard drive) or optical (e.g., a compact disk read only memory, or "CD ROM"), and may be 
read only or random access. Similarly, the transmission medium may be twisted wire pairs, 
coaxial cable, optical fiber, or some other suitable transmission medium known to the art. 
The invention is not limited by these aspects of any given implementation. 
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The particular embodiments disclosed above are illustrative only, as the invention 
may be modified and practiced in different but equivalent manners apparent to those skilled 
in the art having the benefit of the teachings herein. Furthermore, no limitations are intended 
5 to the details of construction or design herein shown, other than as described in the claims 

below. It is therefore evident that the particular embodiments disclosed above may be 
altered or modified and all such variations are considered within the scope and spirit of the 
invention. Accordingly, the protection sought herein is as set forth in the claims below. 
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